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ABSTRACT 

Dissipative axisymmetric pulsar magnetosphere is calculated by a direct numerical simulation of 
the Strong-Field Electrodynamics equations. The magnetic separatrix disappears, it is replaced by a 
region of enhanced dissipation. With a better numerical scheme, one should be able to calculate the 
bolometric lightcurves for a given conductivity. 
Subject headings: pulsars - neutron stars - magnetospheres 



1. INTRODUCTION 

Recent progress in calculating pulsar magnetospheres 
(Spitkovsky 2006 and references therein) is based on 
Force- Free Electrodynamics (FFE). FFE describes elec- 
tromagnetic fields of special geometry - at each event, 
the electric field is smaller than and perpendicular to the 
magnetic field. When applied to pulsars, FFE predicts a 
singular current layer along the magnetic separatrix sur- 
face. Within the light cylinder, the singularity is of an 
admissible type - although the current is singular, the 
electric field remains smaller than and perpendicular to 
the magnetic field. Outside the light cylinder, the singu- 
larity actually violates the FFE equations. For this rea- 
son, and also to understand dissipation processes leading 
to the observed radiation, one wants to use a dissipative 
generalization of FFE. Such a theory - Strong-Field Elec- 
trodynamics (SFE) - has recently been proposed (Gruzi- 
nov 2008). 




Fig. 2. — Left: Invariant Eq, intervals Eq max/5. Right: In- 
variant E 2 — B 2 , intervals (E 2 — B 2 )m ax /5, showing only positive 
isolines. 





Fig. 3. — Left: Poloidal current distribution over the surface of 
the star. Right: Poynting luminosity at different radii. 

SFE describes electromagnetic fields of arbitrary ge- 
ometry. SFE is Maxwell theory with the semi-dissipative 
Lorentz-covariant Ohm's law 



P E x B + (p 2 + 7 2 ct 2 £ 2 ) 1/2 (£qB + £ E) 
B 2 + El 



(1) 



Fig. 1. — a = 10, r B = 0.5. Thick: magnetic stream function 
ip, intervals of i/>(r = 1,0 = 7r/2)/5. Thin: electric potential <j> } 
same intervals. Dotted: poloidal current /, see Fig. 3 for current 
distribution on the surface of the star. 



where p = V-E is the charge density, B^ — Eq = B 2 — E 2 , 
BqEq = E ■ B, Eq > represent the field invariants, and 

7 2 = .} g . The conductivity scalar a = ct(Bq,Eo) is 

an arbitrary function of the field invariants. SFE comes 
from the following simple postulate: at each event, in the 
frame where the electric field is parallel to the magnetic 
field and the charge density vanishes, the current aE 
flows along the fields. *. Covariant formulation of SFE 
is j 2 = -gt 2 £: 2 , B Fj = E Fj (compare to FFE: E = 0. 
Fj=0). 

1 SFE handles time-like 4-currents (Gruzinov 2008). 
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the radius of the star. We then evolved the three com- 
ponents of the electric field (E r ,Eg,E^), the magnetic 
stream function ip, and the toroidal magnetic field B^. 
The boundary conditions on the surface of the star are 
ip = sin 2 9/r s and E e = - sm(26»)/r 2 . 

We show the magnetosphere at t — 6, when the field 
configuration within r = 2 (two light cylinders) satu- 
rates. The outer boundary was at r = 5 — out of causal 
contact with the r < 2 zone. The saturated electric field 
is potential. The saturated toroidal field is given by the 
poloidal current B^ = 21/ (r sin 8). 




Fig. 6.— Same as fig.3, a = 50, r s = 0.2. 



Fig. 4. — Same as fig.l, a = 50, r s 



0.2. 




Fig. 5. — Same as fig.2, cr = 50, r s 



0.2. 



2. AXISYMMETRIC PULSAR IN SFE 

We simulated time-dependent axisymmetric Maxwell 
equations with the SFE Ohm's law. We used the reg- 
ularization described in (Gruzinov 2008), and checked 
that the results are independent of the regularization to 
the specified accuracy. 

The initial field was vacuum around the rotating mag- 
netized neutron star. In spherical coordinates, in pulsar 
units (jj, = il = c = 1), the vacuum field is given by 
the magnetic stream function ip — sin 8/r and the elec- 
trostatic potential <p = (1/3 — cos 2 9) r^/r 3 , where r s is 



We show two cases. The r s = 0.5, cr = 10 results (for 
the Poynting power) are accurate to about 2% . Fig.l. 
shows the field lines, Fig.2. - the field invariants, Fig.3. 
- the integrated poloidal current as a function of the 
latitude and the Poynting luminosity (in pulsar units) as 
a function of radius. The luminosity decreases due to 
dissipation. Figg. 4-6 show the r s = 0.2, cr = 50 results; 
these are accurate to about 10% . 

3. CONCLUSION 

Our code is amateurish - a straightforward discretiza- 
tion with diffusive stabilization 2 . A professional numer- 
ical scheme might be of interest for two reasons: 

(i) Our pulsar, even at cr = 50, is too dissipative. It 
would be interesting to increase a to the values when the 
return current becomes space-like. 

(ii) With a good 3D code, it should be possible to calcu- 
late the pulsar bolometric lightcurves for a given surface 
field and conductivity 3 . To calculate emission, one as- 
sumes efficient heat to light conversion and strong beam- 
ing along the field lines (in the frame where the fields are 
parallel) . 

I thank Andrew MacFadyen for useful advice. This 
work was supported by the David and Lucile Packard 
foundation. 



2 Available upon request from ag92@nyu.edu. 

3 In SFE, cr is an arbitrary function of the field invariants. In 



the real world, the star breaks Lorcntz invariance. 
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